3.3.59 \(\int \frac {1}{x^3 (a+b x^2) (c+d x^2)^3} \, dx\) [259]

3.3.59.1 Optimal result
3.3.59.2 Mathematica [A] (verified)
3.3.59.3 Rubi [A] (verified)
3.3.59.4 Maple [A] (verified)
3.3.59.5 Fricas [B] (verification not implemented)
3.3.59.6 Sympy [F(-1)]
3.3.59.7 Maxima [B] (verification not implemented)
3.3.59.8 Giac [B] (verification not implemented)
3.3.59.9 Mupad [B] (verification not implemented)

3.3.59.1 Optimal result

Integrand size = 22, antiderivative size = 178 \[ \int \frac {1}{x^3 \left (a+b x^2\right ) \left (c+d x^2\right )^3} \, dx=-\frac {1}{2 a c^3 x^2}+\frac {d^2}{4 c^2 (b c-a d) \left (c+d x^2\right )^2}+\frac {d^2 (3 b c-2 a d)}{2 c^3 (b c-a d)^2 \left (c+d x^2\right )}-\frac {(b c+3 a d) \log (x)}{a^2 c^4}+\frac {b^4 \log \left (a+b x^2\right )}{2 a^2 (b c-a d)^3}-\frac {d^2 \left (6 b^2 c^2-8 a b c d+3 a^2 d^2\right ) \log \left (c+d x^2\right )}{2 c^4 (b c-a d)^3} \]

output
-1/2/a/c^3/x^2+1/4*d^2/c^2/(-a*d+b*c)/(d*x^2+c)^2+1/2*d^2*(-2*a*d+3*b*c)/c 
^3/(-a*d+b*c)^2/(d*x^2+c)-(3*a*d+b*c)*ln(x)/a^2/c^4+1/2*b^4*ln(b*x^2+a)/a^ 
2/(-a*d+b*c)^3-1/2*d^2*(3*a^2*d^2-8*a*b*c*d+6*b^2*c^2)*ln(d*x^2+c)/c^4/(-a 
*d+b*c)^3
 
3.3.59.2 Mathematica [A] (verified)

Time = 0.28 (sec) , antiderivative size = 171, normalized size of antiderivative = 0.96 \[ \int \frac {1}{x^3 \left (a+b x^2\right ) \left (c+d x^2\right )^3} \, dx=\frac {1}{4} \left (-\frac {2}{a c^3 x^2}+\frac {d^2}{c^2 (b c-a d) \left (c+d x^2\right )^2}+\frac {2 d^2 (3 b c-2 a d)}{c^3 (b c-a d)^2 \left (c+d x^2\right )}-\frac {4 (b c+3 a d) \log (x)}{a^2 c^4}-\frac {2 b^4 \log \left (a+b x^2\right )}{a^2 (-b c+a d)^3}-\frac {2 d^2 \left (6 b^2 c^2-8 a b c d+3 a^2 d^2\right ) \log \left (c+d x^2\right )}{c^4 (b c-a d)^3}\right ) \]

input
Integrate[1/(x^3*(a + b*x^2)*(c + d*x^2)^3),x]
 
output
(-2/(a*c^3*x^2) + d^2/(c^2*(b*c - a*d)*(c + d*x^2)^2) + (2*d^2*(3*b*c - 2* 
a*d))/(c^3*(b*c - a*d)^2*(c + d*x^2)) - (4*(b*c + 3*a*d)*Log[x])/(a^2*c^4) 
 - (2*b^4*Log[a + b*x^2])/(a^2*(-(b*c) + a*d)^3) - (2*d^2*(6*b^2*c^2 - 8*a 
*b*c*d + 3*a^2*d^2)*Log[c + d*x^2])/(c^4*(b*c - a*d)^3))/4
 
3.3.59.3 Rubi [A] (verified)

Time = 0.40 (sec) , antiderivative size = 174, normalized size of antiderivative = 0.98, number of steps used = 4, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.136, Rules used = {354, 99, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {1}{x^3 \left (a+b x^2\right ) \left (c+d x^2\right )^3} \, dx\)

\(\Big \downarrow \) 354

\(\displaystyle \frac {1}{2} \int \frac {1}{x^4 \left (b x^2+a\right ) \left (d x^2+c\right )^3}dx^2\)

\(\Big \downarrow \) 99

\(\displaystyle \frac {1}{2} \int \left (-\frac {b^5}{a^2 (a d-b c)^3 \left (b x^2+a\right )}-\frac {d^3 \left (6 b^2 c^2-8 a b d c+3 a^2 d^2\right )}{c^4 (b c-a d)^3 \left (d x^2+c\right )}+\frac {-b c-3 a d}{a^2 c^4 x^2}-\frac {d^3 (3 b c-2 a d)}{c^3 (b c-a d)^2 \left (d x^2+c\right )^2}-\frac {d^3}{c^2 (b c-a d) \left (d x^2+c\right )^3}+\frac {1}{a c^3 x^4}\right )dx^2\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {1}{2} \left (\frac {b^4 \log \left (a+b x^2\right )}{a^2 (b c-a d)^3}-\frac {d^2 \left (3 a^2 d^2-8 a b c d+6 b^2 c^2\right ) \log \left (c+d x^2\right )}{c^4 (b c-a d)^3}-\frac {\log \left (x^2\right ) (3 a d+b c)}{a^2 c^4}+\frac {d^2 (3 b c-2 a d)}{c^3 \left (c+d x^2\right ) (b c-a d)^2}+\frac {d^2}{2 c^2 \left (c+d x^2\right )^2 (b c-a d)}-\frac {1}{a c^3 x^2}\right )\)

input
Int[1/(x^3*(a + b*x^2)*(c + d*x^2)^3),x]
 
output
(-(1/(a*c^3*x^2)) + d^2/(2*c^2*(b*c - a*d)*(c + d*x^2)^2) + (d^2*(3*b*c - 
2*a*d))/(c^3*(b*c - a*d)^2*(c + d*x^2)) - ((b*c + 3*a*d)*Log[x^2])/(a^2*c^ 
4) + (b^4*Log[a + b*x^2])/(a^2*(b*c - a*d)^3) - (d^2*(6*b^2*c^2 - 8*a*b*c* 
d + 3*a^2*d^2)*Log[c + d*x^2])/(c^4*(b*c - a*d)^3))/2
 

3.3.59.3.1 Defintions of rubi rules used

rule 99
Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_)*((e_.) + (f_.)*(x_) 
)^(p_), x_] :> Int[ExpandIntegrand[(a + b*x)^m*(c + d*x)^n*(e + f*x)^p, x], 
 x] /; FreeQ[{a, b, c, d, e, f, p}, x] && IntegersQ[m, n] && (IntegerQ[p] | 
| (GtQ[m, 0] && GeQ[n, -1]))
 

rule 354
Int[(x_)^(m_.)*((a_) + (b_.)*(x_)^2)^(p_.)*((c_) + (d_.)*(x_)^2)^(q_.), x_S 
ymbol] :> Simp[1/2   Subst[Int[x^((m - 1)/2)*(a + b*x)^p*(c + d*x)^q, x], x 
, x^2], x] /; FreeQ[{a, b, c, d, p, q}, x] && NeQ[b*c - a*d, 0] && IntegerQ 
[(m - 1)/2]
 

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 
3.3.59.4 Maple [A] (verified)

Time = 2.78 (sec) , antiderivative size = 187, normalized size of antiderivative = 1.05

method result size
default \(-\frac {1}{2 a \,c^{3} x^{2}}+\frac {\left (-3 a d -b c \right ) \ln \left (x \right )}{a^{2} c^{4}}-\frac {b^{4} \ln \left (b \,x^{2}+a \right )}{2 a^{2} \left (a d -b c \right )^{3}}+\frac {d^{3} \left (-\frac {c^{2} \left (a^{2} d^{2}-2 a b c d +b^{2} c^{2}\right )}{2 d \left (d \,x^{2}+c \right )^{2}}+\frac {\left (3 a^{2} d^{2}-8 a b c d +6 b^{2} c^{2}\right ) \ln \left (d \,x^{2}+c \right )}{d}-\frac {c \left (2 a^{2} d^{2}-5 a b c d +3 b^{2} c^{2}\right )}{d \left (d \,x^{2}+c \right )}\right )}{2 c^{4} \left (a d -b c \right )^{3}}\) \(187\)
norman \(\frac {-\frac {1}{2 a c}+\frac {\left (6 a^{2} d^{3}-10 a b c \,d^{2}+3 b^{2} c^{2} d \right ) d \,x^{4}}{2 a \,c^{3} \left (a^{2} d^{2}-2 a b c d +b^{2} c^{2}\right )}+\frac {\left (9 a^{2} d^{3}-15 a b c \,d^{2}+4 b^{2} c^{2} d \right ) d^{2} x^{6}}{4 c^{4} a \left (a^{2} d^{2}-2 a b c d +b^{2} c^{2}\right )}}{\left (d \,x^{2}+c \right )^{2} x^{2}}-\frac {b^{4} \ln \left (b \,x^{2}+a \right )}{2 a^{2} \left (a^{3} d^{3}-3 a^{2} b c \,d^{2}+3 a \,b^{2} c^{2} d -b^{3} c^{3}\right )}-\frac {\left (3 a d +b c \right ) \ln \left (x \right )}{a^{2} c^{4}}+\frac {d^{2} \left (3 a^{2} d^{2}-8 a b c d +6 b^{2} c^{2}\right ) \ln \left (d \,x^{2}+c \right )}{2 c^{4} \left (a^{3} d^{3}-3 a^{2} b c \,d^{2}+3 a \,b^{2} c^{2} d -b^{3} c^{3}\right )}\) \(297\)
risch \(\frac {-\frac {d^{2} \left (3 a^{2} d^{2}-5 a b c d +b^{2} c^{2}\right ) x^{4}}{2 c^{3} a \left (a^{2} d^{2}-2 a b c d +b^{2} c^{2}\right )}-\frac {d \left (9 a^{2} d^{2}-15 a b c d +4 b^{2} c^{2}\right ) x^{2}}{4 c^{2} a \left (a^{2} d^{2}-2 a b c d +b^{2} c^{2}\right )}-\frac {1}{2 a c}}{\left (d \,x^{2}+c \right )^{2} x^{2}}-\frac {3 \ln \left (x \right ) d}{a \,c^{4}}-\frac {\ln \left (x \right ) b}{a^{2} c^{3}}-\frac {b^{4} \ln \left (b \,x^{2}+a \right )}{2 a^{2} \left (a^{3} d^{3}-3 a^{2} b c \,d^{2}+3 a \,b^{2} c^{2} d -b^{3} c^{3}\right )}+\frac {3 d^{4} \ln \left (-d \,x^{2}-c \right ) a^{2}}{2 c^{4} \left (a^{3} d^{3}-3 a^{2} b c \,d^{2}+3 a \,b^{2} c^{2} d -b^{3} c^{3}\right )}-\frac {4 d^{3} \ln \left (-d \,x^{2}-c \right ) a b}{c^{3} \left (a^{3} d^{3}-3 a^{2} b c \,d^{2}+3 a \,b^{2} c^{2} d -b^{3} c^{3}\right )}+\frac {3 d^{2} \ln \left (-d \,x^{2}-c \right ) b^{2}}{c^{2} \left (a^{3} d^{3}-3 a^{2} b c \,d^{2}+3 a \,b^{2} c^{2} d -b^{3} c^{3}\right )}\) \(396\)
parallelrisch \(-\frac {-9 x^{6} a^{4} d^{6}-6 a^{3} b \,c^{4} d^{2}+6 a^{2} b^{2} c^{5} d -32 \ln \left (x \right ) x^{2} a^{3} b \,c^{3} d^{3}+24 \ln \left (x \right ) x^{2} a^{2} b^{2} c^{4} d^{2}+16 \ln \left (d \,x^{2}+c \right ) x^{2} a^{3} b \,c^{3} d^{3}-32 \ln \left (x \right ) x^{6} a^{3} b c \,d^{5}+32 \ln \left (d \,x^{2}+c \right ) x^{4} a^{3} b \,c^{2} d^{4}-12 \ln \left (d \,x^{2}+c \right ) x^{2} a^{2} b^{2} c^{4} d^{2}+24 \ln \left (x \right ) x^{6} a^{2} b^{2} c^{2} d^{4}+16 \ln \left (d \,x^{2}+c \right ) x^{6} a^{3} b c \,d^{5}-12 \ln \left (d \,x^{2}+c \right ) x^{6} a^{2} b^{2} c^{2} d^{4}-64 \ln \left (x \right ) x^{4} a^{3} b \,c^{2} d^{4}+48 \ln \left (x \right ) x^{4} a^{2} b^{2} c^{3} d^{3}-24 \ln \left (d \,x^{2}+c \right ) x^{4} a^{2} b^{2} c^{3} d^{3}-2 a \,b^{3} c^{6}+2 a^{4} c^{3} d^{3}+24 x^{6} a^{3} b c \,d^{5}-19 x^{6} a^{2} b^{2} c^{2} d^{4}+4 x^{6} a \,b^{3} c^{3} d^{3}+32 x^{4} a^{3} b \,c^{2} d^{4}-26 x^{4} a^{2} b^{2} c^{3} d^{3}+6 x^{4} a \,b^{3} c^{4} d^{2}-4 \ln \left (x \right ) x^{6} b^{4} c^{4} d^{2}+2 \ln \left (b \,x^{2}+a \right ) x^{6} b^{4} c^{4} d^{2}+24 \ln \left (x \right ) x^{4} a^{4} c \,d^{5}-8 \ln \left (x \right ) x^{4} b^{4} c^{5} d +4 \ln \left (b \,x^{2}+a \right ) x^{4} b^{4} c^{5} d -12 \ln \left (d \,x^{2}+c \right ) x^{4} a^{4} c \,d^{5}+12 \ln \left (x \right ) x^{2} a^{4} c^{2} d^{4}-6 \ln \left (d \,x^{2}+c \right ) x^{2} a^{4} c^{2} d^{4}-12 x^{4} a^{4} c \,d^{5}+12 \ln \left (x \right ) x^{6} a^{4} d^{6}-6 \ln \left (d \,x^{2}+c \right ) x^{6} a^{4} d^{6}-4 \ln \left (x \right ) x^{2} b^{4} c^{6}+2 \ln \left (b \,x^{2}+a \right ) x^{2} b^{4} c^{6}}{4 \left (a^{3} d^{3}-3 a^{2} b c \,d^{2}+3 a \,b^{2} c^{2} d -b^{3} c^{3}\right ) \left (d \,x^{2}+c \right )^{2} x^{2} a^{2} c^{4}}\) \(675\)

input
int(1/x^3/(b*x^2+a)/(d*x^2+c)^3,x,method=_RETURNVERBOSE)
 
output
-1/2/a/c^3/x^2+(-3*a*d-b*c)/a^2/c^4*ln(x)-1/2*b^4/a^2/(a*d-b*c)^3*ln(b*x^2 
+a)+1/2*d^3/c^4/(a*d-b*c)^3*(-1/2*c^2*(a^2*d^2-2*a*b*c*d+b^2*c^2)/d/(d*x^2 
+c)^2+(3*a^2*d^2-8*a*b*c*d+6*b^2*c^2)/d*ln(d*x^2+c)-c*(2*a^2*d^2-5*a*b*c*d 
+3*b^2*c^2)/d/(d*x^2+c))
 
3.3.59.5 Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 640 vs. \(2 (168) = 336\).

Time = 5.04 (sec) , antiderivative size = 640, normalized size of antiderivative = 3.60 \[ \int \frac {1}{x^3 \left (a+b x^2\right ) \left (c+d x^2\right )^3} \, dx=-\frac {2 \, a b^{3} c^{6} - 6 \, a^{2} b^{2} c^{5} d + 6 \, a^{3} b c^{4} d^{2} - 2 \, a^{4} c^{3} d^{3} + 2 \, {\left (a b^{3} c^{4} d^{2} - 6 \, a^{2} b^{2} c^{3} d^{3} + 8 \, a^{3} b c^{2} d^{4} - 3 \, a^{4} c d^{5}\right )} x^{4} + {\left (4 \, a b^{3} c^{5} d - 19 \, a^{2} b^{2} c^{4} d^{2} + 24 \, a^{3} b c^{3} d^{3} - 9 \, a^{4} c^{2} d^{4}\right )} x^{2} - 2 \, {\left (b^{4} c^{4} d^{2} x^{6} + 2 \, b^{4} c^{5} d x^{4} + b^{4} c^{6} x^{2}\right )} \log \left (b x^{2} + a\right ) + 2 \, {\left ({\left (6 \, a^{2} b^{2} c^{2} d^{4} - 8 \, a^{3} b c d^{5} + 3 \, a^{4} d^{6}\right )} x^{6} + 2 \, {\left (6 \, a^{2} b^{2} c^{3} d^{3} - 8 \, a^{3} b c^{2} d^{4} + 3 \, a^{4} c d^{5}\right )} x^{4} + {\left (6 \, a^{2} b^{2} c^{4} d^{2} - 8 \, a^{3} b c^{3} d^{3} + 3 \, a^{4} c^{2} d^{4}\right )} x^{2}\right )} \log \left (d x^{2} + c\right ) + 4 \, {\left ({\left (b^{4} c^{4} d^{2} - 6 \, a^{2} b^{2} c^{2} d^{4} + 8 \, a^{3} b c d^{5} - 3 \, a^{4} d^{6}\right )} x^{6} + 2 \, {\left (b^{4} c^{5} d - 6 \, a^{2} b^{2} c^{3} d^{3} + 8 \, a^{3} b c^{2} d^{4} - 3 \, a^{4} c d^{5}\right )} x^{4} + {\left (b^{4} c^{6} - 6 \, a^{2} b^{2} c^{4} d^{2} + 8 \, a^{3} b c^{3} d^{3} - 3 \, a^{4} c^{2} d^{4}\right )} x^{2}\right )} \log \left (x\right )}{4 \, {\left ({\left (a^{2} b^{3} c^{7} d^{2} - 3 \, a^{3} b^{2} c^{6} d^{3} + 3 \, a^{4} b c^{5} d^{4} - a^{5} c^{4} d^{5}\right )} x^{6} + 2 \, {\left (a^{2} b^{3} c^{8} d - 3 \, a^{3} b^{2} c^{7} d^{2} + 3 \, a^{4} b c^{6} d^{3} - a^{5} c^{5} d^{4}\right )} x^{4} + {\left (a^{2} b^{3} c^{9} - 3 \, a^{3} b^{2} c^{8} d + 3 \, a^{4} b c^{7} d^{2} - a^{5} c^{6} d^{3}\right )} x^{2}\right )}} \]

input
integrate(1/x^3/(b*x^2+a)/(d*x^2+c)^3,x, algorithm="fricas")
 
output
-1/4*(2*a*b^3*c^6 - 6*a^2*b^2*c^5*d + 6*a^3*b*c^4*d^2 - 2*a^4*c^3*d^3 + 2* 
(a*b^3*c^4*d^2 - 6*a^2*b^2*c^3*d^3 + 8*a^3*b*c^2*d^4 - 3*a^4*c*d^5)*x^4 + 
(4*a*b^3*c^5*d - 19*a^2*b^2*c^4*d^2 + 24*a^3*b*c^3*d^3 - 9*a^4*c^2*d^4)*x^ 
2 - 2*(b^4*c^4*d^2*x^6 + 2*b^4*c^5*d*x^4 + b^4*c^6*x^2)*log(b*x^2 + a) + 2 
*((6*a^2*b^2*c^2*d^4 - 8*a^3*b*c*d^5 + 3*a^4*d^6)*x^6 + 2*(6*a^2*b^2*c^3*d 
^3 - 8*a^3*b*c^2*d^4 + 3*a^4*c*d^5)*x^4 + (6*a^2*b^2*c^4*d^2 - 8*a^3*b*c^3 
*d^3 + 3*a^4*c^2*d^4)*x^2)*log(d*x^2 + c) + 4*((b^4*c^4*d^2 - 6*a^2*b^2*c^ 
2*d^4 + 8*a^3*b*c*d^5 - 3*a^4*d^6)*x^6 + 2*(b^4*c^5*d - 6*a^2*b^2*c^3*d^3 
+ 8*a^3*b*c^2*d^4 - 3*a^4*c*d^5)*x^4 + (b^4*c^6 - 6*a^2*b^2*c^4*d^2 + 8*a^ 
3*b*c^3*d^3 - 3*a^4*c^2*d^4)*x^2)*log(x))/((a^2*b^3*c^7*d^2 - 3*a^3*b^2*c^ 
6*d^3 + 3*a^4*b*c^5*d^4 - a^5*c^4*d^5)*x^6 + 2*(a^2*b^3*c^8*d - 3*a^3*b^2* 
c^7*d^2 + 3*a^4*b*c^6*d^3 - a^5*c^5*d^4)*x^4 + (a^2*b^3*c^9 - 3*a^3*b^2*c^ 
8*d + 3*a^4*b*c^7*d^2 - a^5*c^6*d^3)*x^2)
 
3.3.59.6 Sympy [F(-1)]

Timed out. \[ \int \frac {1}{x^3 \left (a+b x^2\right ) \left (c+d x^2\right )^3} \, dx=\text {Timed out} \]

input
integrate(1/x**3/(b*x**2+a)/(d*x**2+c)**3,x)
 
output
Timed out
 
3.3.59.7 Maxima [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 364 vs. \(2 (168) = 336\).

Time = 0.23 (sec) , antiderivative size = 364, normalized size of antiderivative = 2.04 \[ \int \frac {1}{x^3 \left (a+b x^2\right ) \left (c+d x^2\right )^3} \, dx=\frac {b^{4} \log \left (b x^{2} + a\right )}{2 \, {\left (a^{2} b^{3} c^{3} - 3 \, a^{3} b^{2} c^{2} d + 3 \, a^{4} b c d^{2} - a^{5} d^{3}\right )}} - \frac {{\left (6 \, b^{2} c^{2} d^{2} - 8 \, a b c d^{3} + 3 \, a^{2} d^{4}\right )} \log \left (d x^{2} + c\right )}{2 \, {\left (b^{3} c^{7} - 3 \, a b^{2} c^{6} d + 3 \, a^{2} b c^{5} d^{2} - a^{3} c^{4} d^{3}\right )}} - \frac {2 \, b^{2} c^{4} - 4 \, a b c^{3} d + 2 \, a^{2} c^{2} d^{2} + 2 \, {\left (b^{2} c^{2} d^{2} - 5 \, a b c d^{3} + 3 \, a^{2} d^{4}\right )} x^{4} + {\left (4 \, b^{2} c^{3} d - 15 \, a b c^{2} d^{2} + 9 \, a^{2} c d^{3}\right )} x^{2}}{4 \, {\left ({\left (a b^{2} c^{5} d^{2} - 2 \, a^{2} b c^{4} d^{3} + a^{3} c^{3} d^{4}\right )} x^{6} + 2 \, {\left (a b^{2} c^{6} d - 2 \, a^{2} b c^{5} d^{2} + a^{3} c^{4} d^{3}\right )} x^{4} + {\left (a b^{2} c^{7} - 2 \, a^{2} b c^{6} d + a^{3} c^{5} d^{2}\right )} x^{2}\right )}} - \frac {{\left (b c + 3 \, a d\right )} \log \left (x^{2}\right )}{2 \, a^{2} c^{4}} \]

input
integrate(1/x^3/(b*x^2+a)/(d*x^2+c)^3,x, algorithm="maxima")
 
output
1/2*b^4*log(b*x^2 + a)/(a^2*b^3*c^3 - 3*a^3*b^2*c^2*d + 3*a^4*b*c*d^2 - a^ 
5*d^3) - 1/2*(6*b^2*c^2*d^2 - 8*a*b*c*d^3 + 3*a^2*d^4)*log(d*x^2 + c)/(b^3 
*c^7 - 3*a*b^2*c^6*d + 3*a^2*b*c^5*d^2 - a^3*c^4*d^3) - 1/4*(2*b^2*c^4 - 4 
*a*b*c^3*d + 2*a^2*c^2*d^2 + 2*(b^2*c^2*d^2 - 5*a*b*c*d^3 + 3*a^2*d^4)*x^4 
 + (4*b^2*c^3*d - 15*a*b*c^2*d^2 + 9*a^2*c*d^3)*x^2)/((a*b^2*c^5*d^2 - 2*a 
^2*b*c^4*d^3 + a^3*c^3*d^4)*x^6 + 2*(a*b^2*c^6*d - 2*a^2*b*c^5*d^2 + a^3*c 
^4*d^3)*x^4 + (a*b^2*c^7 - 2*a^2*b*c^6*d + a^3*c^5*d^2)*x^2) - 1/2*(b*c + 
3*a*d)*log(x^2)/(a^2*c^4)
 
3.3.59.8 Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 357 vs. \(2 (168) = 336\).

Time = 0.29 (sec) , antiderivative size = 357, normalized size of antiderivative = 2.01 \[ \int \frac {1}{x^3 \left (a+b x^2\right ) \left (c+d x^2\right )^3} \, dx=\frac {b^{5} \log \left ({\left | b x^{2} + a \right |}\right )}{2 \, {\left (a^{2} b^{4} c^{3} - 3 \, a^{3} b^{3} c^{2} d + 3 \, a^{4} b^{2} c d^{2} - a^{5} b d^{3}\right )}} - \frac {{\left (6 \, b^{2} c^{2} d^{3} - 8 \, a b c d^{4} + 3 \, a^{2} d^{5}\right )} \log \left ({\left | d x^{2} + c \right |}\right )}{2 \, {\left (b^{3} c^{7} d - 3 \, a b^{2} c^{6} d^{2} + 3 \, a^{2} b c^{5} d^{3} - a^{3} c^{4} d^{4}\right )}} + \frac {18 \, b^{2} c^{2} d^{4} x^{4} - 24 \, a b c d^{5} x^{4} + 9 \, a^{2} d^{6} x^{4} + 42 \, b^{2} c^{3} d^{3} x^{2} - 58 \, a b c^{2} d^{4} x^{2} + 22 \, a^{2} c d^{5} x^{2} + 25 \, b^{2} c^{4} d^{2} - 36 \, a b c^{3} d^{3} + 14 \, a^{2} c^{2} d^{4}}{4 \, {\left (b^{3} c^{7} - 3 \, a b^{2} c^{6} d + 3 \, a^{2} b c^{5} d^{2} - a^{3} c^{4} d^{3}\right )} {\left (d x^{2} + c\right )}^{2}} - \frac {{\left (b c + 3 \, a d\right )} \log \left (x^{2}\right )}{2 \, a^{2} c^{4}} + \frac {b c x^{2} + 3 \, a d x^{2} - a c}{2 \, a^{2} c^{4} x^{2}} \]

input
integrate(1/x^3/(b*x^2+a)/(d*x^2+c)^3,x, algorithm="giac")
 
output
1/2*b^5*log(abs(b*x^2 + a))/(a^2*b^4*c^3 - 3*a^3*b^3*c^2*d + 3*a^4*b^2*c*d 
^2 - a^5*b*d^3) - 1/2*(6*b^2*c^2*d^3 - 8*a*b*c*d^4 + 3*a^2*d^5)*log(abs(d* 
x^2 + c))/(b^3*c^7*d - 3*a*b^2*c^6*d^2 + 3*a^2*b*c^5*d^3 - a^3*c^4*d^4) + 
1/4*(18*b^2*c^2*d^4*x^4 - 24*a*b*c*d^5*x^4 + 9*a^2*d^6*x^4 + 42*b^2*c^3*d^ 
3*x^2 - 58*a*b*c^2*d^4*x^2 + 22*a^2*c*d^5*x^2 + 25*b^2*c^4*d^2 - 36*a*b*c^ 
3*d^3 + 14*a^2*c^2*d^4)/((b^3*c^7 - 3*a*b^2*c^6*d + 3*a^2*b*c^5*d^2 - a^3* 
c^4*d^3)*(d*x^2 + c)^2) - 1/2*(b*c + 3*a*d)*log(x^2)/(a^2*c^4) + 1/2*(b*c* 
x^2 + 3*a*d*x^2 - a*c)/(a^2*c^4*x^2)
 
3.3.59.9 Mupad [B] (verification not implemented)

Time = 6.68 (sec) , antiderivative size = 314, normalized size of antiderivative = 1.76 \[ \int \frac {1}{x^3 \left (a+b x^2\right ) \left (c+d x^2\right )^3} \, dx=-\frac {\frac {1}{2\,a\,c}+\frac {x^4\,\left (3\,a^2\,d^4-5\,a\,b\,c\,d^3+b^2\,c^2\,d^2\right )}{2\,a\,c^3\,\left (a^2\,d^2-2\,a\,b\,c\,d+b^2\,c^2\right )}+\frac {x^2\,\left (9\,a^2\,d^3-15\,a\,b\,c\,d^2+4\,b^2\,c^2\,d\right )}{4\,a\,c^2\,\left (a^2\,d^2-2\,a\,b\,c\,d+b^2\,c^2\right )}}{c^2\,x^2+2\,c\,d\,x^4+d^2\,x^6}-\frac {\ln \left (d\,x^2+c\right )\,\left (3\,a^2\,d^4-8\,a\,b\,c\,d^3+6\,b^2\,c^2\,d^2\right )}{-2\,a^3\,c^4\,d^3+6\,a^2\,b\,c^5\,d^2-6\,a\,b^2\,c^6\,d+2\,b^3\,c^7}-\frac {b^4\,\ln \left (b\,x^2+a\right )}{2\,\left (a^5\,d^3-3\,a^4\,b\,c\,d^2+3\,a^3\,b^2\,c^2\,d-a^2\,b^3\,c^3\right )}-\frac {\ln \left (x\right )\,\left (3\,a\,d+b\,c\right )}{a^2\,c^4} \]

input
int(1/(x^3*(a + b*x^2)*(c + d*x^2)^3),x)
 
output
- (1/(2*a*c) + (x^4*(3*a^2*d^4 + b^2*c^2*d^2 - 5*a*b*c*d^3))/(2*a*c^3*(a^2 
*d^2 + b^2*c^2 - 2*a*b*c*d)) + (x^2*(9*a^2*d^3 + 4*b^2*c^2*d - 15*a*b*c*d^ 
2))/(4*a*c^2*(a^2*d^2 + b^2*c^2 - 2*a*b*c*d)))/(c^2*x^2 + d^2*x^6 + 2*c*d* 
x^4) - (log(c + d*x^2)*(3*a^2*d^4 + 6*b^2*c^2*d^2 - 8*a*b*c*d^3))/(2*b^3*c 
^7 - 2*a^3*c^4*d^3 + 6*a^2*b*c^5*d^2 - 6*a*b^2*c^6*d) - (b^4*log(a + b*x^2 
))/(2*(a^5*d^3 - a^2*b^3*c^3 + 3*a^3*b^2*c^2*d - 3*a^4*b*c*d^2)) - (log(x) 
*(3*a*d + b*c))/(a^2*c^4)